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The nonequilibrium dynamic phase transition, in the ki- 
netic Ising model in presence of an oscillating magnetic field, 
has been studied by Monte Carlo simulation. The fluctuation 
of dynamic order parameter has been studied as a function of 
temperature near the dynamic transition point. The tempera- 
ture variation of appropriately defined 'susceptibility' has also 
been studied near the dynamic transition point. Similarly, 
fluctuation of energy and appropriately defined 'specific-heat' 
have been studied as a function of temperature near the dy- 
namic transition point. In both the cases, the fluctuations (of 
dynamic order parameter and energy) and the corresponding 
responses diverge (as power law fashion) near the dynamic 
transition point with similar critical behavior (with identical 
exponent values). 



PACS number(s): 05.50.+q 



I. INTRODUCTION 

The physics of equilibrium phase transition in the Ising 
model is well understood jjj. However, the mechanism 
behind the nonequilibrium phase transition is not yet ex- 
plored rigorously and the basic phenomenology is still 
undeveloped. It is quite interesting to study how the 
system behaves if it is driven out of equilibrium. The 
simplest prototype example is the kinetic Ising model in 
oscillating magnetic field. In this context, the dynamic 
response of the Ising system in presence of an oscillat- 
ing magnetic field has been studied extensively &-|7j in 
the last few years. The dynamic hysteresis 0411 and 
the nonequilibrium dynamic phase transition Iq-R] are 
two important aspects of the dynamic response of the 
kinetic Ising model in presence of an oscillating mag- 
netic field. The nonequilibrium dynamic phase transition 
in the kinetic Ising model in presence of an oscillating 
magnetic field, was first studied by Tome and Oliviera 
H. They solved the mean- field (MF) dynamic equation 
of motion (for the average magnetization) of the kinetic 
Ising model in presence of a sinusoidally oscillating mag- 
netic field. By defining the dynamic order parameter as 
the time averaged magnetization over a full cycle of the 
oscillating magnetic field, they showed that depending 



upon the value of the field amplitude and the tempera- 
ture, the dynamic order parameter takes nonzero value 
from a zero value. In the field amplitude and temperature 
plane there exists a distinct phase boundary separating 
dynamic ordered (nonzero value of order parameter) and 
disordered (order parameter vanishes) phase. A tricriti- 
cal point (TCP), (separating the nature (discontinuous- 
continuous) of the transition)on the phase boundary line, 
was also observed by them However, one may argue 
that such a mean-field transition is not truly dynamic 
in origin since it exists even in the quasi-static (or zero 
frequency) limit. This is because, if the field amplitude 
is less than the coercive field (at temperature less than 
the transition temperature without any field), then the 
response magnetization varies periodically but asymmet- 
rically even in the zero frequency limit; the system re- 
mains locked to one well of the free energy and cannot 
go to the other one, in the absence of fluctuation. 

The true dynamic nature of this kind of phase tran- 
sition (in presence of fluctuation) was first attemted to 
study by Lo and Pelcovits Q. They have studied the 
dynamic phase transition in the kinetic Ising model in 
presence of an oscillating magnetic field by Monte Carlo 
(MC) simulation which allows the microscopic fluctua- 
tions. Here, the transition disappears in the zero fre- 
quency limit; due to the fluctuations, the magnetization 
flips to the direction of the magnetic field and the dy- 
namic order parameter (time averaged magnetization) 
vanishes. However, they fa] have not reported any precise 
phase boundary. Acharyya and Chakrabarti studied 
the nonequilibrium dynamic phase transition in the ki- 
netic Ising model in presence of oscillating magnetic field 
by extensive MC simulation. They have also iden- 
tified that this dynamic phase transition (at a particu- 
lar nonzero frequency of the oscillating magnetic field) 
is associated with the breaking of the symmetry of the 
dynamic hysteresis (m — h) loop. In the dynamically 
disordered phase (where the value of order parameter 
vanishes) the corresponding hysteresis loop is symmet- 
ric, and loses its symmetry in the ordered phase (givin 
nonzero value of dynamic order parameter). They 
also studied the temperature variation of the ac suscep- 
tibility components near the dynamic transition point. 
The major observation was that the imaginary (real) part 
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of the ac susceptibility gives a peak (dip) near the dy- 
namic transition point (where the dynamic order param- 
eter vanishes). The important conclusions were: (i) this 
is a distinct signal of phase transition and (ii) this is an in- 
dication of the thermodynamic nature of the phase tran- 
sition. The Debye relaxation of the dynamic order pa- 
rameter and the critical slowing down have been studied 
very recently |]llf both by MC simulation and by solving 
the dynamic MF equation || of motion for the average 
magnetization. The specific-heat singularity near the 
dynamic transition point is also an indication of the ther- 
modynamic nature of this dynamic phase transition. It 
is worth mentioning here that the statistical distribution 
of dynamic order parameter has been studied by Sides 
et al ||. The nature of the distribution changes (from 
bimodal to unimodal) near the dynamic transition point. 
They have also observed || that the fluctuation of the 
hysteresis loop area becomes considerably large near the 
dynamic transition point. 

In the case of equilibrium phase transitions, the fluc- 
tuation - dissipation theorem (FDT) states (due to the 
applicability of Gibbs formalism) that the mean square 
fluctuations of some intrinsic physical quantities (say, en- 
ergy, magnetization etc.) are directly related with some 
responses (specific heat, susceptibility etc.) of the sys- 
tem. Consequently, near the ferro-para transition point, 
both the fluctuation of magnetization and the suscepti- 
bilty show same singular behavior. If it is of power law 
type, the same singular behavior will be characterised by 
the same exponent. This is also true for fluctuation of en- 
ergy and the specific heat. These are the consequences of 
fluctuation-dissipation theorem [0. Here, the main mo- 
tivation is to study the fluctuations and corresponding 
responses near the dynamic transition temperature. 

In this paper, the fluctuations of dynamic order pa- 
rameter and the energy are studied as a function of tem- 
perature near the dynamic transition point. The temper- 
ature variations of 'susceptibility' and the 'specific- heat' 
are also studied near the transition point. The temper- 
ature variation of the fluctuation of dynamic order pa- 
rameter and that of the 'susceptibility' are compared. 
Similarly, the temperature variation of the fluctuation of 
energy and that of the 'specific-heat' are compared. The 
paper is organised as follows: the model and the simu- 
lation scheme are discussed in section II, the results are 
reported in Section III, section IV contains the summary 
of the work. 



II. MODEL AND SIMULATION 

The Ising model with nearest neighbor ferromagnetic 
coupling in presence of a time varying magnetic field can 
be represented by the Hamiltonian 

H = ~ £ Jij^-h®^ (2.1) 

<ij> i 



Here, sf(= ±1) is Ising spin variable, Jy is the in- 
teraction strength and h(t) — hocos(u)t) represents the 
oscillating magnetic field, where ho and lo are the am- 
plitude and the frequency respectively of the oscillating 
field. The system is in contact with an isothermal heat 
bath at temperature T . For simplicity all Jy (= J > 0) 
are taken equal to unity and the boundary condition is 
chosen to be periodic. The temperature (T) is measured 
in the unit of J/K, where K is the Boltzmann constant 
(here K is taken unity). 

A square lattice of linear size L(= 100) has been con- 
sidered. At any finite temperature T and for a fixed 
frequency (u>) and amplitude (ho) of the field, the mi- 
croscopic dynamics of this system has been studied here 
by Monte Carlo simulation using Glauber single spin-flip 
dynamics with a particular choice of the Metropolis rate 
of single spin- flip fl^| . Starting from an initial condition 
where all spins are up, each lattice site is updated here 
sequentially and one such full scan over the entire lat- 
tice is defined as the unit time step (Monte Carlo step 
or MCS). The instanteneous magnetization (per site), 
m(t) — (l/L 2 )J2i s i nas been calculated. From the in- 
stanteneous magnetization, the dynamic order parame- 
ter Q = £ § m(t)dt (time averaged magnetization over 
a full cycle of the oscillating field) is calculated. Some of 
the transient loops have been discarded to get the stable 
value of the dynamical quantities. 

III. RESULTS 

A. Temperature variations of susceptibility and 
fluctuation of dynamic order parameter 

The fluctuation of the dynamic order parameter is 

SQ 2 = (< Q 2 > - < Q > 2 ) , 

where the <> stands for the averaging over various 
Monte Carlo samples. 

The 'susceptibility' is defined as 

_ d <Q > 
^ dho 

Here, a square lattice of linear size L (=100) has been 
considered. < Q 2 > and < Q > are calculated using 
MC simulation. The averaging has been done over 100 
different (uncorrelated) MC samples. 

The temperature variations of fluctuation of Q, i.e., 
5Q 2 and 'susceptibility' \ have been studied here and 
both plotted in Fig. 1. From the figure it is observed 
that both SQ 2 and x diverge near the dynamic transition 
point (where Q vanishes). 

This has been studied for two different values of field 
amplitude ho (Fig. la is for ho = 0.2 and Fig. lb is for 
ho =0.1). The dynamic transition temperatures Td(ho), 
at which x and SQ 2 diverge, are 1.91±0.01 for ho = 0.2 
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and 2.15 ± 0.01 for h = 0.1. These values of T d (h ) 
agree with the phase diagram estimated from vanishing 
of Q. The log e (x) versus log e (Td — T) and log e (SQ 2 ) 
versus log e (T d — T) plots show (insets of Fig. 1) that 
X ~ (T d - T)- a and 6Q 2 ~ (T d - T)-<*. For h = 0.2, 
a ~ 0.53 (inset of Fig. la) and for ho = 0.1, a ~ 2.5 
(inset of Fig. lb). Results show that both \ an d 5Q 2 
diverge near T d as a power law with the same exponent 
a, though there is a crossover region (where the effective 
exponent values are different). 



B. Temperature variations of specific- heat and 
fluctuation of energy 

The time averaged (over a full cycle) cooperative en- 
ergy of the system is 

E=-(u;/2nL 2 )<f ^l>>',j 

and the fluctuation of the cooperative energy is 

SE 2 = (< E 2 > - < E > 2 ) . 

The 'specific-heat' C |ll[] is defined as the derivative of 
the energy (defined above) with respect to the tempera- 
ture, and 



C 



d< E> 
dT ■ 



Here, also a square lattice of linear size L (=100) has 
been considered. < E 2 > and < E > are calculated 
using MC simulation. The averaging has been done over 
100 different (uncorrelated) MC samples. 

The temperature variation of the 'specific-heat' has 
been studied fll|| and found prominent divergent behav- 
ior near the dynamic transition point (where < Q > van- 
ishes). The temperature variation of fluctuation of en- 
ergy, (8E) 2 has been studied and plotted in Fig. 2. From 
the figure it is clear that the mean square fluctuation of 
energy (SE 2 ) and the specific heat (C) both diverge near 
the dynamic transition point (where the dynamic order 
parameter Q vanishes) 

This has been studied for two different values of field 
amplitude ho (Fig. 2a is for ho = 0.2 and Fig. 2b is for ho 
=0.1). Here also (like the earlier case) the specific heat 
C and SE 2 are observed to diverge at temperatures T d . 
The temperatures T d (ho), at which C and 8E 2 diverge, 
are 1.91±0.01 for h = 0.2 (Fig. 2a) and 2.15±0.01 for h 
— 0.1 (Fig. 2b). These values also agree with the phase 
diagram estimated from vanishing of Q. The log e (C) 
vs. log e (T d - T) and log e {SE 2 ) vs. log e (T d - T) plots 
show (insets of Fig. 2) that C ~ (T d - T)~ 7 and 5E 2 ~ 
(T d -T)-T. For h a = 0.2, 7 - 0.35 (inset of Fig. 2a) and 
for ho = 0.1, 7 ~ 0.43 (inset of Fig. 2b). Like the earlier 



case, here also the results show that both \ and SQ 2 
diverge near T d as a power law with the same exponent 
a, though there is a crossover region (where the effective 
exponent values are different). 



IV. SUMMARY 

The nonequlibrium dynamic phase transition, in the 
kinetic Ising model in presence of an oscillating magnetic 
field, is studied by Monte Carlo simulation. Acharyya 
and Chakrabarti @ observed that the complex suscepti- 
bility components have peaks (or dips) at the dynamic 
transition point. Sides et al || observed that the fluctu- 
ation in the hysteresis loop area grows (seems to diverge) 
near the dynamic transition point. It has been observed 
|TTf that the 'relaxation time' and the appropriately de- 
fined 'specific-heat' diverge near the dynamic transition 
point. 

The mean square fluctuation of dynamic order param- 
eter and the 'susceptibility' are studied as a function of 
temperature, near the dynamic transition point. Both 
shows the power law variation with respect to the re- 
duced temperature near the dynamic transition point 
with the same exponent values. Similar, observation has 
been made for the case of mean square fluctuation of en- 
ergy and the 'specific-heat'. It appears that although the 
effective exponent values for the fluctuation and the ap- 
propriate linear resoponse differ considerably, away from 
the dynamic transition point T d , they eventually con- 
verge and give identical value as the temperature inter- 
val \T d — T\ decreases and falls within a narrow crossover 
region. These numerical observations indicate that the 
fluctuation- dissipation relation jjj holds good in this case 
of the nonequilibrium phase transition in the kinetic Ising 
model. However, at this stage there is no analytic sup- 
port of the FDT in this case. 

Finally, it should be mentioned, in this context, that 
experiments [|| on ultrathin ferromagnetic Fe/Au(001) 
films have been performed to study the frequency depen- 
dence of hysteresis loop areas. Recently, attempts have 
been made pX[ | to measure the dynamic order parame- 
ter Q experimentally, in the same material, by extend- 
ing their previous study |J. The dynamic phase tran- 
sition has been studied from the observed temperature 
variation of Q. However, the detailed investigation of 
the dynamic phase transitions by measuring variations 
of associated response functions (like the ac susceptibil- 
ity, specific-heat, correlations, relaxations etc) have not 
yet been performed experimentally. 
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(triangles). All the data points are plotted in arbitrary 
units. Here, u> = 2nx 0.01. Corresponding insets show 
the plots of log e (C) (circles) and log e (SE 2 ) (triangles) 
against log e (Td~T). Solid lines represent the linear best- 
fit in the region very closed to X^. 
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Figure Captions 



Fig.l Temperature variations of < Q >, x an d 5Q 2 for 
two different values of field amplitudes ho'- (a) For h = 
0.2; Q (solid line), \ (circles) and SQ 2 (triangles), (b) 
For ho = 0.1; Q (solid line), \ (circles) and 5Q 2 (trian- 
gles). All the data points are plotted in arbitrary units. 
Here, u> = 2irx 0.01. Corresponding insets show the 
plots of log e (x) (circles) and log e (SQ 2 ) (triangles) against 
log e (Td — T). Solid lines represent the linear best-fit in 
the region very closed to T4. 



Fig. 2 Temperature variations of < Q >, C and 5E 2 for 
two different values of field amplitudes ho: (a) For ho 
= 0.2; Q (solid line), C (circles) and 5E 2 (triangles), 
(b) For h = 0.1; Q (solid line), C (circles) and SE 2 
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